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Previous studies indicate that nonlinear properties of Gaussian time series with long-range corre- 
lations, Ui, can be detected and quantified by studying the correlations in the magnitude series \ui\, 
i.e., the "volatility". However, the origin for this empirical observation still remains unclear, and 
the exact relation between the correlations in Ui and the correlations in \ui\ is still unknown. Here 
we find analytical relations between the scaling exponent of linear series Ui and its magnitude series 
\ui\. Moreover, we find that nonlinear time series exhibit stronger (or the same) correlations in the 
magnitude time series compared to linear time series with the same two-point correlations. Based 
on these results we propose a simple model that generates multifractal time series by explicitly 
inserting long range correlations in the magnitude series; the nonlinear multifractal time series is 
generated by multiplying a long-range correlated time series (that represents the magnitude series) 
with uncorrelated time series [that represents the sign series sgn(ui)\. Our results of magnitude 
series correlations may help to identify linear and nonlinear processes in experimental records. 

PACS numbers: 87.10.-l-e, 89.20.-a, 89.65.Gh, 89.75.Da 
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I. INTRODUCTION 

Natural systems often exhibit irregular and complex 
behavior that at first look erratic but in fact posses scale 
invariant structure [e.g.,[3,|2|- In many cases this nontriv- 
ial structure points to long-range temporal correlations, 
which means that very far events are actually (statisti- 
cally) correlated with each other. Long-range correla- 
tions are usually characterized by scaling laws where the 
scaling exponents quantify the strength of these correla- 
tions. However, it is clear that the two-point long-range 
correlations reveal just one aspect of the complexity of 
the system under consideration and that higher order 
statistics is needed to fully characterized the statistical 
properties of the system. 

The two-point correlation function is usually used to 
quantify the scale invariant structure of time series (long- 
range correlations), while the g-point correlation func- 
tion quantifies the higher order correlations. In some 
cases the g-point correlation function is trivially related 
to the two-point correlation function, and the scaling ex- 
ponents of different moments are linearly dependent on 
the second moment scaling exponent. These kind of pro- 
cesses are termed "linear" and "mono fractal" , since just 
a single exponent that determines the two-point corre- 
lations (and thus the linear correlations) quantifies the 
entire spectrum of q order scaling exponents. In other 
cases, the relation between the g-point correlation func- 
tion has nontrivial relation to the two-point correlation 
function, and a (nontrivial) spectrum of scaling expo- 
nents is needed to quantify the statistical properties of 
the system; these processes are called "nonlinear" and 
"multifractal" . The classification into linear and nonlin- 
ear processes is important for the understandings of the 
underlying dynamics of natural time series and for mod- 
els development. Moreover, the nonlinear properties of 
natural time series may have practical diagnosis use [e.g.. 



Recently a simple measure for nonlinearity of time se- 
ries was suggested Q- Given a time series Ui, the cor- 
relations in the magnitude series (volatility) |ui| may be 
related (in some cases) to the nonlinear properties of the 
time series; basically, when the magnitude series is cor- 
related the time series Ui is nonlinear. It was also shown 
that the scaling exponent of the magnitude series may 
be related, in some cases, to the multifractal spectrum 
width. However, these observations are empirical and the 
reasons underly these observations still remain unclear. 

Here we develop an analytical relation between the 
scaling exponent of the original time series Ui and the 
scaling exponent of the magnitude time series \ui\ for 
linear series. We show that when the original time se- 
ries is nonlinear, the corresponding scaling exponent of 
the magnitude series is larger (or in some cases equal) 
compared to that of linear series and that the correla- 
tions in the magnitude series increase as the nonlinear- 
ity of the original series increases. These relations may 
help to identify nonlinear processes and quantify the non- 
linearity strength. Based on these results we suggest a 
generic model for multifractality by multiplying random 
signs with long-range correlated noise, and show that the 
multifractal spectrum width and the volatility exponent 
increase as these correlations become stronger. 

The paper is organized as follows: in Section ^ we 
present some background regarding non-linear processes 
and magnitude (volatility) series correlations. In Section 
mil we develop an analytical relation between the orig- 
inal time series scaling exponent a and the magnitude 
series exponent a^] we confirm the analytical relation 
using numerical simulation. We then study in Section lTvl 
the relation between volatility correlations and the mul- 
tifractal spectrum width of several multifractal models, 
and introduce a simple model that generates multifrac- 
tal series by explicitly inserting long range correlations in 



the magnitude series. A summary of the results is g: 
in Section [V] 
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II. NONLINEARITY AND VOLATILITY 
CORRELATIONS 

A. Two-point correlations 

The long range correlations of a time series {ui} 
{i = 0,1,2, ... ,N) can be evaluated using the two-point 
correlation function (uiUj) ((•) stands for expectation 
value); when Ui is long-range correlated and stationary 
the two-point correlation function is (uiUj) ~ |i — j\~^'^ 
(0 < 7 < 1) jj, y|. It is possible to get good estima- 
tion of the scaling exponent using various methods, such 
as the power spectrum, Fluctuation Analysis (FA) [g, 
Detrended Fluctuation Analysis (DFA) @, 0, Q , wavelet 
transform Q, and others; see |j| for more details. The 
different techniques characterize the linear two point cor- 
relations in a time series with a scaling exponent which 
is related to the scaling exponent 7. 

In this study we use the FA method for the analyt- 
ical derivations since this method is relatively simple. 
In the FA method the sequence is treated as steps of 
a random walk; then the variance of its displacement, 
Xt = X]i=o^»' ^^ found by averaging over different time 
windows of length t. The scaling exponent a of the series 
(also referred to as the Hurst exponent H) can be mea- 
sured using the relation Var{Xt) = {Xt^) ^ t^" where 
Var{-) is the variance; the scaling exponent a is related 
to the correlation exponent 7 by 2 — 7 = 2a. 



B. High order correlations 

A more complete description of stochastic process {ui} 
with a zero mean is given by its multivariate distribution: 
P{uq, Ui, U2, ...). It is equivalent to the knowledge stored 
in the correlation functions of different orders |3, Hojl : 
{ui), (uiUj), (uiUjUk), (uiUjUkUi), etc. In many cases 
it is useful to use the cummulants of different orders Cq 
which are related to the q order correlation function by 
113: 



Ci 


= M 


(1) 


C2 


= (UiUj) 


(2) 


Cs 


= (UiUjUk) 


(3) 


Ci 


= (UiUjUkUl) - {UiUj)(ukUl) - 






{UiUk){UjUl) ~ {u^Ul){UjUk), 


(4) 



and so on. 

For a linear process (sometimes referred to as "Gaus- 
sian" process), all cummulants above the second are 
equal to zero (Wicks theorem) ^^. Thus, in this 
case, the two-point correlation fully describes the pro- 
cess [III, u^, since all correlation functions (of positive 



and even order) may be expressed as products of the 
two-point correlation function (uiUj). 

Processes that are nonlinear (or "multifractal" ) have 
nonzero high order cummulants. The nonlinearity of 
these processes may be detected by measuring the multi- 
fractal spectrum jUO using advanced techniques, such 
as the Wavelet Transform Modulus Maxima p| or the 
Multifractal DFA (MF-DFA) _14|. In MF-DFA we calcu- 
late the q order correlation function of the profile Xt — 
^j^P Ui and the partition function is Zq{t) = {\Xt\'^). 
For time series that obey scaling laws the partition func- 
tion is Zq{t) - t9"(9). Thus, the "spectrum" of scahng 
exponents a(q) characterizes the correlation functions of 
different orders. For a linear series, the exponents a{q) 
will all give a single value a for all q [lj|. 



C. Volatility correlations 

A known example for the use of volatility correlations 
(defined below) is econometric time series Il9|. Econo- 
metric time series exhibit irregular behavior such that the 
changes (logarithmic increments) in the time series have 
a white noise spectrum (uncorrelated). Nonetheless, the 
magnitudes of the changes exhibit long-range correlations 
that reflect the fact that economic markets experience 
quiet periods with clusters of less pronounced price fluc- 
tuations (up and down), followed by more volatile periods 
with pronounced fluctuations (up and down). This type 
of correlation is referred to as "volatility correlations" . 

Given a time series Ui, the magnitude (volatility) se- 
ries may be defined as \A.Ui\ = |wi+i — Ui\. The scaling 
exponent of the magnitude series is the volatility scaling 
exponent ay. Gorrelations in the magnitude series are 
observed to be closely related to nonlinearity and multi- 
fractality pi, 16, 17] 

In this paper we refer to "volatility" with two small 
differences. First, we consider the square of the series 
elements rather than their absolute values. According to 
our abservations, this transformation has negligible effect 
on the scaling exponent ay, but it substantially simplifies 
the analytical treatment. Second, for simplicity, we also 
consider the series itself rather than the increment series. 
That is: the volatility series is defined as uf rather than 
lAwij. Note that in most applications the absolute values 
of the increment series are considered instead of the ab- 
solute values of the series itself, since the original series is 
mostly nonstationary (defined below); here we overcome 
this problem by first considering stationary series. 



D. Stationary and nonstationary time series 

Series with < a < 1 are stationary, that is, their cor- 
relation function depends only on the difference between 
points i and j, i.e., (uiUj) = f{\i — j\). Their variance is a 
finite constant, that does not increase with the sequence 
length. Sequences with a > 1 are non stationary and 
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FIG. 1: Magnitude series scaling exponent Qi, vs. the two- 
point correlation exponent a for linear sequences Ui. The solid 
line represents results for synthesized sequences of length 2^^, 
averaged over 15 configurations, for uf and \ui\ (these two 
coincide). The circles represent the analytical reconstruction 
taking into account corrections due to finite size effects and 
non-stationarity. Analytical results for N -^ oo are given by 
the dashed line. 



have a different form of correlation function where the 
correlation function depends also on the absolute indices 
i and j, {u^UJ) = z2«-2+j2"-2_|i _ j|2"-2. ggg ^^ gcal- 

ing exponents of nonstationary series (or series with poly- 
nomial trends) may be calculated using methods that can 
eliminate constant or polynomial trends from the data 



III. VOLATILITY CORRELATIONS OF LINEAR 
TIME SERIES 

We proceed to study the relation between the volatility 
correlation exponent a^ and the original scaling exponent 
a for linear processes, both numerically and analytically. 



A. Simulations 

We generate artificial long-range correlated linear se- 
quences Ui with different values of a in the range a G 
(0, 1.5] as follows 18]: (i) generate Gaussian white noise 
series, (ii) apply Fourier transform on that series, (iii) 
multiply the power spectrum S{f) by 1/ f^ where (3 = 
2a — 1 and / > 0, and (iv) apply inverse Fourier trans- 
form. The resultant series is long-range correlated with 
a scaling exponent a. We measure the volatility scaling 
exponent a^, i.e. the scaling exponent oiuf and \ui\, ver- 
sus the original scaling exponent. The results are plotted 
in Fig. n 

The simulations indicate that the dependence of a^ on 



a for linear series may be divided into three regions: for 
a < 3/4 we obtain ay « 1/2, for a > 1.25 we obtain 
ay ~ a, while for 0.75 < a < 1.25 there is a transi- 
tion region between those two behaviors. These results 
were obtained using the DFA method which can handle 
nonstationary time series. 

It is important to note that in Fig. ^ we used the ab- 
solute value of the series Ui, to calculate the volatility 
scaling exponent. Nevertheless, the usual method for cal- 
culating the volatility scaling exponent is performed by 
taking the absolute value of the differences series, Am^, 
rather than Ui itself [3, [l3 • The reason for this is that in 
many cases the given time series has an exponent in the 
range ^ < a < l^. By differentiating the sequence we 
get a new sequence Aui with an exponent a = a — 1 < ^ . 
According to our analysis, if the sequence is linear, the 
volatility exponent for this series will be a^, = ^ (whereas 
for the original series a^ may be higher than i even for 
linear data). Thus, if the scaling exponent of the magni- 
tude of the differences series, \Aui\, differs from 1/2 this 
is an indication for nonlinearity. 

We note that for a ^ 1 the series is highly non- 
stationary, i.e., it is most of the time either above or 
below 0, apart from a discrete set of crossing points. The 
behavior of the series Ui is not very different than the 
behavior of its absolute value \ui\, and therefore it is not 
surprising that a^ = a. 



B. Analytical Treatment 

Let us consider a Gaussian distributed linear sequence 
Ui of length t with scaling exponent a. For simplicity, 
we assume that the sequence is stationary. Consider the 
magnitude series: Vi = uf. In order to calculate the 
magnitude series scaling exponent a„ we will calculate 
the variance of the displacement Vt — X]i=o ^*- 



Var{Vt) = {Vt^)-{Vt){Vt) 



ELoELoKf^vj) - {'"t){vj)] 



= ELoEj=o[K'"j> 



K'>K)]- 



Because the series Ui is linear, the fourth cumulant is 
C4 — (Wicks theorem), and by using Eq. Q) we get, 



{u^u'^ = {u^){u'^+2{u,u,)^ (5) 



and thus, 



i=o j=a 

Substituting the two-point correlation function for long- 
range correlated data: 

p{^^J) = {u,u,)^!^|-^^~' ;^^; (6) 



we obtain 



-27 



Var{Vt)^t + J2\i-j\ 

^t + t-27+2. 

Since 2 — 7 = 2a the above expression becomes 



(F,2) ^t + t" 



j-2q„ 



(7) 



(8) 



For a < I the first term, i, is dominant and for i — > 00 
we obtain a„ = i. Otherwise the second term, i**"^^, is 
dominating and thus a^ ~ 2a — 1. However, the sim- 
ulation results (Fig. ^ deviate from at, w 2a — 1 as 
a ^ 1. This is because as a — *■ 1, logarithmic and poly- 
nomial corrections due to strong finite size effects and 
non-stationarity must be taken into account in our cal- 
culations (i.e., the variance of the sequence depends on 
its length; see Appendix ^ . This is done by dividing 
{V^) [Eq. ©] by the variance of the sequence: 



Var{ui) = 



1 



1-a 



^(l-i2"-2) 



const 
Int 



t 



2a-2 



a< 1 

a = 1 (9) 
a> 1 



This modification yields an a^ that is very close to the 
one obtained from the numerical simulation (in the tran- 
sition region 0.75 < a < 1.25, and also for a > 1.25 with 
tti, — a; see Fig. ^. The relation a^ = a for a > 1.25 
can now be proved analytically: It is noticeable that the 
dominant scaling term of (V^) for the nonstationary case 
is proportional to i'*""^ j^gq^ jj^jj Dividing by the vari- 
ance term, t^"~^ [Eq. Q], yields t'^" ^ t^"" and hence 
a — a,,- 



IV. VOLATILITY CORRELATIONS AND THE 
MULTIFRACTAL SPECTRUM WIDTH 

A. Random multifractal cascades 

Following I3, [135 we study the relation between the 
volatility scaling exponent a„ and the multifractal spec- 
trum width of nonlinear multifractal time series. We gen- 
erate artificial noise with multifractal properties accord- 
ing to the algorithm proposed in [16j. The algorithm is 
based on random cascades on wavelet dyadic trees. The 
multifractal series is constructed by building its wavelet 
coefficients at different scales recursively, where at each 
stage the coefficients of the coarser scale are multiplied 
by a random variable W in order to built the coefficients 
of the finer scale. Note that we now consider the in- 
crements series, hence the generated time series is sta- 
tionary. The multifractal spectrum /(a) depends on the 
statistical properties of the random variable W. 

We choose W to follow log-normal distribution, such 
that the ln|W^| is normally distributed with fj, and a^ 



being the mean and variance. For this case the multi- 
fractal spectrum /(a) is known analytically |16J| and by 
assigning /(a) = it possible to obtain a„ 



'^rmn^raax -> 



V2cr 



a„ 



In 2 
In 2' 



(10) 
(11) 



Thus, 

2 v^g- 



the multifractal width, Aa = amax — ctmin = 
, depends just on a while the scaling exponent a(0) 
depends on fi, i.e. a(0) — —-^ [iq . 

Using the above algorithm, we generate multifractal 
time series with a fixed multifractal width Aa (by fixing 
a) and different scaling exponents a{q — 2) (by chang- 
ing /i), and calculate their volatility exponents ay (see 
Fig. O [13 ■ We perform the same analysis for their 
respective surrogate time series, which are linearized se- 
ries (after phase randomization) that have the same two- 
point correlations with exponent a{q = 2) as the original 
series |20|; note that a{q) is the g-point correlation ex- 
ponent. We find that the volatility exponent a^, calcu- 
lated in section nm for the linear case, is the lower bound 
for all multifractal sequences (studied here) with same 
a{q = 2). Nonetheless, for a{q = 2) > 1, i.e., for non- 
stationary series, a^ = a{q = 2) as in linear series. It 
is clearly seen in Fig. |3 that for stationary time series 
(a(2) < 1), the volatility correlations increase as the mul- 
tifractal spectrum width becomes wider, or alternatively, 
as the nonlinearity of the original series strengthens. 



B. A simple model for multifractality 

We now propose a simple model for generating multi- 
fractal records, based on the property that multifractal 
series exhibit long range correlations in the volatility se- 
ries. Following '2l'|, we multiply a long range correlated 
series rji (with a scaling exponent a,, > 0.75) with a se- 
ries of uncorrelated random signs e^ = ±1. The resultant 
series, Ui = eirji, has a two-point correlation exponent 
a = 1/2 because of the random signs e^. The magnitude 
exponent a„ is the same as the magnitude exponent ay^ri 
for r]i, because \ui\ = \r]i\. Thus, using our results from 
Section ITTTl if we take a^ > 0.75 we get a sequence with 
a — 1/2 and a^ ~ 2a,, — 1 > 1/2 (see Fig.|21 full diamond 
symbol for a,, — 0.95). Note that in Fig. |2]the theoret- 
ical value of a,, = 2 X 0.95 — 1 = 0.9 is higher than that 
of the numerical estimation a„ — 0.8, most probably due 
to finite size effects. 

According to our derivation in Section IIIII this se- 
quence is nonlinear/ multifractal (because linear series 
with a two-point correlation exponent a — 1/2 should 
have ttv — 1/2). Indeed, one can see from Fig. O that 
the multifractal width for this model increases as a,, in- 
creases beyond 0.75. 

Natural processes are often characterized by complex 
nonlinear and multifractal properties. However, the un- 
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FIG. 2: The magnitude series exponent a,, vs. the two-point 
correlation exponent a for multifractal series. The full tri- 
angles and squares represent sequences generated by the log- 
normal random cascade algorithm with a = 0.1 (triangles) 
and a = 0.05 (squares). The respective linear (phase ran- 
domized) surrogate data sequences are represented by empty 
symbols. The solid line indicates simulation results for lin- 
ear sequences as derived in [13, explained in section ITTTl and 
shown in Fig. The full diamond represents the scaling ex- 
ponent of our multifractal model's sequences Ui — et-qi with 
an = 0.95, while the empty diamond represents the respec- 
tive scaling exponent of the surrogate data. All sequences 
are of length 2^^ elements, and results were averaged over 15 
configurations. Error bars are smaller than symbol size. 



FIG. 3; (a) Multifractal spectrum width and (b) volatility 
exponent a„ for sequences of the form Ui — eirji of length 
2^^, averaged over 15 configurations. The error bars indicate 
the mean ± 1 std. For a^ > 0.75 both volatility correlation 
exponent and the multifractal spectrum width of the series 
are increasing with a,,. 



resent the relatively slow oceanic processes. Our model 
can also be used to model other complex systems like 
economy and heartbeat dynamics. 



SUMMARY 



derlying mechanisms of these processes are usually not 
so well understood. Several prototypes for multifractal 
processes include, e.g., (i) the energy cascade model de- 
scribing turbulence [a Il2| , (ii) the universal multifractal 
process usually used to generally explain geophysical phe- 
nomena 22, 2^], and (iii) the turbulence-like model for 
heart rate variability |24| . 

The multifractal model described in this section is a 
simple model with known properties that may help to 
gain better understanding of multifractal processes. The 
model consists of two components as follows, (i) a random 
series (which can be also any other long-range correlated 
series) that may represent fast processes of a natural sys- 
tem, which as a first approximation may be regarded as a 
white noise, interacting with (ii) a long-range correlated 
process that may represents a slow modulation of the 
natural system. This interaction results in episodes with 
less volatile fluctuations followed by episodes with more 
volatile fluctuations. In the context of heart-rate vari- 
ability, the fast component may represents the parasym- 
pathetic branch of the autonomic nervous system while 
the slow process may represents the sympathetic and the 
hormonal activities. In the context of geophysical phe- 
nomena, the fast cornponent may represents the fast at- 
mospheric processes |25j while the slow process may rep- 



We study the behavior of the magnitude series scaling 
exponent ay versus the original two-point scaling expo- 
nent a for linear and nonlinear (multifractal) series. We 
find analytically and by simulations that for linear series 
the dependence of a„ on a may be divided into three 
regions: for a < 3/4 the volatility exponent is ay = 1/2, 
for a > 1.25 the volatility exponent is Uy = a, while 
for 0.75 < a < 1.25 there is a transition region in which 
logarithmic corrections due to finite size effects and non- 
stationarity are dominant. 

The results presented here provide the theory for the 
relation found previously 0, [13 between multifractality 
and the scaling exponent of the magnitude of the differ- 
ences series (volatility). This relation provides a simple 
method for preliminary detection and quantification of 
nonlinear time series, a procedure which usually requires 
relatively complex techniques. 

We also study the volatility of some known models of 
multifractal time series, and find that their magnitude 
scaling exponent is bounded from below by ay of the 
corresponding phase randomized linear surrogate series. 

Based on the above findings, we propose a simple 
model that generates multifractality by explicitly insert- 
ing long range correlations (a,, > 0.75) into the magni- 
tude series. This model may serve as a generic model for 



multifractality and may help to gain preliminary under- 
standings of natural complex phenomena. The model, 
which involves interaction between fast and slow compo- 
nents, may represents natural fast processes that interact 
with slower processes. In addition, the simplicity of the 
model may help to identify these processes more easily 
in experimental records. 
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APPENDIX A: FINITE SIZE EFFECTS AND 
NON-STATIONARITY NEAR a = 1 

A linear time sequence with scaling exponent a can be 
generated by filtering Gaussian white noise such that the 
power spectrum will be _18J : 



S{f) 



/ = 



(Al) 



where (3 = 2a — 1. Assume a signal Ui of N discrete 
points sampled at time intervals At. The power spectrum 
consists of N points in the frequency range (— 2At' 2Sa\ 
with intervals of A/ = jf^- Thus, looking only at the 
positive frequencies, the minimal frequency (without loss 



of generality) is -if 



The variance of the signal 



2 2NAf 

is the total area under the power spectrum: 



/23T fTSt 1 

Sif)df = 2 / —^df (A2) 
■JWKt "^ 2JvAt ■' 

Assuming At = 1, for a — 1 the variance is, 



Var{u^) = 2 In TV. 



(A3) 



Thus, the variance diverges logarithmically for a = 1. 
For a 7^ 1 the variance is 



Var{ui) = 



1 



■)2a-2 



1 



(1 - iv^"-^). 



(A4) 



From ljA4|) follows: For a < 1 the variance converges, 
and for a > 1 it diverges. 

Non-stationarity: For a > 1 the variance diverges 
with the sequence length N, because of the singularity in 
the power spectrum, and the sequence is non-stationary. 
For a > 1 the divergence is power-law, i.e. Var(ui) ^ 
7V2q-2^ while at a = 1 the divergence is logarithmic. 

Finite size effects: For a < 1 the variance converges 
to a finite constant so the sequence is stationary, but as 



a — > 1 this convergence becomes slower. This means 
that as a — > 1, larger and larger sequence lengths N 
are required so that the variance will indeed converge 
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FIG. 4: Correlation coefficient p(0) (i.e. the variance {u^}) for 
linear sequences of A^ = 50, 000 points. Circles indicate the 
simulation results. Dots represent analytical results for the 
variance calculated according to Eq. (IA4II . which takes into 
account the finite size effects. The solid line is the variance 
for N —^ oo. It can be seen that as a -^ 1 the convergence 
becomes slower and finite size effects become more dominant 
[i.e., the convergence is non-uniform in the range a £ (0, 1)]. 



to a constant value (see Fig. 0J. This argument also 
holds for other values of the correlation functions p{n), 
n = 0, 1, ...,oo, although in a more moderate way. 

The strong finite size effects around 0.75 < a < 1.25 
and the non-stationarity at a > 1 have to be taken into 
account when calculating the magnitude series scaling 
exponent a^ . This is done by dividing the volatility fluc- 
tuation function (V^) by the variance of the sequence 
given in Eq. HA4|I [2g | . 

For N —f oo the finite size effects disappear and ay 
converges to its theoretical value (see Fig.QJ. This con- 
vergence is extremely slow and becomes weaker as we 
approach a ~ 1. 

For completeness, we show in Fig. \^ the correlation 
coefficients p{n = i — j) for a < 1/2. In this regime the 
sequences exhibits short range anti-correlations as can be 
seen in Fig. [S] The expression of the correlation function 
for a < 1/2 is approximately Q 



p{i - j) = (u^Uj) 



a{2a 
2a 



m-j\ 



2a-2 



i ¥" j 



(A5) 
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FIG. 5: Correlation coefficients p(0) (circles) and p(l) 
(squares) for linear sequences of 50,000 points, in range 
a < 1/2. The dashed line indicates the analytical results for 
p(0), taking into account the finite series size effects, which 
approximately follows results for N ^ ca (solid line). p(l) is 
negative for a < 1/2 indicating anti-correlations. The solid 
lines are the analytical expressions of Eq. lASt . 
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